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PROPERTIES  OF  EVEN-LENGTH  BARKER  CODES 
AND  SPECIFIC  POLYPHASE  CODES  WITH 
BARKER  TYPE  AUTOCORRELATION  FUNCTION 


INTRODUCTION 

A  Barker  code  is  a  sequence  of  A  numbers  x,  (where  x,  *■>  ±1),  which  has  the  following  auto¬ 
correlation  function: 

N_K  A  for  X  -  0 

/UX)«T  X,x,+k-  (la) 

0  or  ±  1  for  X  -  1 ,  2 . (A  -  1) 

i.e..  the  "time  sidelobes"  in  the  autocorrelation  function  do  not  exceed  the  level  of  1. 

In  radar  applications,  the  sequence  modulates  the  phase  of  the  signal  (some  constant  carrier  fre¬ 
quency)  from  code  element  to  code  element.  For  a  stationary  target  the  above  property  (la)  holds,  but 
now,  since  the  time  variable  is  continuous,  we  pet  small  triangles  in  the  autocorrelation  function  whose 
peaks  are  0  or  ±  1,  and  a  big  triangle  whose  peak  is  A  (the  match  point).  For  a  moving  target  we  actu¬ 
ally  have  the  crosscorrclation  function  of  the  transmitted  code  and  the  target  return,  resulting  in  higher 
sidelobes.  Only  the  autocorrelation  function  will  be  considered  here. 

The  known  code  lengths  having  the  property  of  Eq.  (la)  are  2,  3,  4,  5,  7,  1 1,  13  (I). 

It  has  been  shown  that  no  Barker  code  of  odd  length  exists  for  A'  >  13.  Also,  if  an  even-length 
Barker  code  exists,  it  must  be  a  perfect  square  |2),  i.e.,  A  -  I*.  Since  A  is  even,  /is  also  even. 

The  purpose  here  is  to  investigate  the  possibility  of  even-length  Barker  codes  greater  than  the 
known  of  length  2  (+  +  and  -  +)  and  4  (+■*--  4-  and  +  +  +  -).  Possible  candidates  for  this  are,  for 
example,  lengths  of  16,  36,  64,  100,  etc.,  but  is  was  verified  (2J  that  up  to  A**  6084  (/ —  78)  no 
Barker  code  exists. 

If  x,  is  not  restricted  to  +1,  -I,  but  can  be  any  complex  number  whose  magnitude  Is  unity 
U,|  -  I,  then  the  autocorrelation  function  is  required  to  fulfill: 

A  for  X  -  0 
w-x 

Jj  *♦**♦*“  o  or  ^  unity  magnitude  *  ^ 

for  K  -  1,  2 . A  -  l 

tn  general,  R  (X)  is  a  complex  number.  The  complex  conjugate  is  denoted  by  *. 

DEFINITION  (for  convenience):  A  code  with  property  (lb)  is  a  polphode.  It  is  actually  a 
polyphase  code  with  Barker  type  autocorrelation  function  (excluding  the  real  Barker  codes).  Specific 
types  of  polphodcs  arc  the  generalized  Barker  codes  131  which  are  derived  from  a  "father*  real  Barker 
code.  These  will  be  discussed  later. 

MuuwnjX  tubmiital  Fctxuwy  9,  I  SHI 
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The  following  analysis  will  investigate  the  properties  of  even-length  Barker  codes  and  polphodes 
(where  N  =  I2,  Wand  / are  even),  if  they  exist.  The  analysis  of  Barker  codes  (for  which  Turyn  [2]  con¬ 
siders  evidence  overwhelming  that  they  do  not  exist)  will  lead  to  the  analysis  of  the  general  case  of  pol¬ 
phodes. 

GENERAL  ANALYSIS:  SPECTRUM 

The  general  description  of  a  phase-coded  signal  is  shown  in  Fig.  1.  We  are  interested  in  a  con¬ 
stant  amplitude  code;  thus,  without  loss  of  generality,  we  assume  its  amplitude  is  1,  and  its  carrier  fre¬ 
quency  is  constant  /Q. 


The  signal  duration  7*  is  divided  into  N  code  elements,  each  of  T/N  duration,  and  each  code  ele¬ 
ment  has  phase  <t>,  (for  Barker  codes  <!>,  can  take  only  G.rr  values  corresponding  to  real  x,  which  equal 

+1,  ~ l  in  the  sequence),  where  /  «•  I.  2 . N.  For  polphodes,  <t>,  can  take  any  value  resulting  in  a 

complex  sequence  jr#.  We  will  specify  the  restrictions  on  whenever  they  apply. 

Taking  out  the  carrier  frequency,  the  complex  envelope  of  each  code  element  is  lZ&  “  **''  The 
analysis  from  now  on  will  be  carried  out  with  the  complex  envelope. 


The  spectrum  of  the  signal  is 

sw  -  ij  w>. 

where  $,(/)  is  the  spectrum  of  the  /***  code  clement: 

$</)-/  -  J  * 

m  ■'Xif-n 

$</> «.  Jo  SiMr'***  +  j  r  st(/)r'**0  +  ... 

15 

+  f  f  ,  SN(th~'iul,dt. 

r~i 

and  after  a  change  of  variables  in  the  integrals  (in  order  to  have  the  same  limits  in  each  one) 
Sift  -  /0*  eJ*'c-»’t'di  +  /'**'*//  /V'*** 

+  ^  f*  e,*ie~}3*,>dt  +  .... 


(2) 


0) 


(4) 


(5) 
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S(J)  -  -jLj  eJ*'  (l  -  e-J^rr/N j  +  _  e-j2n/T/N^e~j2n/T/N 

+  _  e-j2vJT/N^e-j2vf2T/N  +  t 

U  +  e«*r*frm 

-Jiirf  \ 


+  ei(*i^nTiN)  +  _ 

tiv  f  1  rl  sin  (2ff/T/ 


U  \  2  N\  hr f-  T/2N  [Terms 


Above 


define 


ln/2 H-* 


and  i-  is  a  scaled  frequency  variable.  Then, 


S{fl  -  -=t|  +  /♦’-* 

2A  4> 

+  /*r*>  +  +  U2*1  . 


This  is  the  basic  spectrum  expression  that  we  will  utilize  through  the  analysis.  The  sin  tfi/tji  term 
in  Eq.  (SO)  is  due  the  basic  code  element  length  77 A.  and  the  terms  in  the  right  bracket  arc  due  to  the 
phase  coding  inside  the  code. 

If  the  signal  bandwidth  is  &  and  we  sample  it  at  lire  Nyquist  rate,  then  T/N  «  1  /B  (this  is 
because  in  general  we  use  /  and  Q  processing,  which  requires  sampling  at  once,  and  not  twice,  the 
reciprocal  of  the  bandwidth).  In  this  case  ii>  **>  vf/B  and  (-*1/2)  T/N  •  -1/2&  But  we  will  proceed 
with  the  general  analysis. 


The  power  spectrum  is 


(£(/)! 2  -  S(J)  S*(A 


and  it  is  the  Fourier  transform  of  the  autocorrelation  function.  Note  that  !$(/}|3  is  always  a  real  func¬ 
tion  of  f,  and  H (r)  is  an  even  function  of  r  for  real  codes,  while  Hit)  *•  /(*(-r)  for  complex  codes. 


To  see  this  relation  in  the  discrete  phase  code,  let  us  examine  in  detail  Barker  codes  of  lengths  7 


and  4. 


BARKER  CODE  7 

This  code  is  known  to  be: 

+  +  +  --  +  - 
0  0  0  tr  w  0  it' 
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Notice  here  that  changing  the  signs  of  all  the  code  elements  does  not  change  the  property  of  the 
autocorrelation  function.  This  means  that  one  can  choose  arbitrarily  the  sign  of  the  first  code  element. 
This  is  true  for  any  Barker  code,  and  polphode,  and  we  will  choose  the  first  code  element  to  be  a'j  =*  4-1 
(or  equivalently  0i  =  0°)  from  now  on,  unless  otherwise  stated. 

T 

The  autocorrelation  function  is  shown  in  Fig.  2  (where  r  =  K-,  K  =  0,  1 . N  -  1). 


t 


MATCH  POINT 
(PEAK  OF  CORRELATION! 

ftf  2  -  Autocorrelation  function  of  Barker  code  7 


According  to  liq.  (10),  substituting  the  known  <t>,  for  this  code  we  get: 


S(/)- 


_L 

2N 


sinjfc. 


j  4.  f-t**  +  4-  t  (12a) 


and 


e#  HL&.  |1+fit*  +  gi**  _  em  _  eiu  4.  gjm  _ 

ti» 


(12b) 


Carrying  out  the  multiplication  of  liqs.  (12a)  and  (12b),  we  get: 

U(/)P  -  si/)  s*(f)  -  ^ 


7  -  t~t**  gj**  _  f-jt*  - 


-—r  ^  7  —  2  cos  4^1  —  2  cos  8gr  -  2  cos  1 2gr  . 

4(VJ  0 


(13) 


We  see  in  Eq.  (13)  the  Fourier  transform  relation  between  |5(/)P  and  the 
lion;  the  7  term  in  the  square  bracket  of  Eq.  (13)  gives  the  ^correlation  peak. 


width  is  one  code  clement  T/N,  is  the  result  of  the 


sin  0 

* 


term,  os  known 


autocorrelation  func- 
(Thc  triangle,  whose 

by  Fourier  transform 
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theory.)  The  three  sidelobes  (on  each  side  of  the  match  point)  are  the  result  of  the  2  cos  4<|<,  2  cos  8t/», 
2  cos  12 0  terms  in  Eq.  (13),  (for  convenience,  we  will  call  these  terms  in  |S(/)P  "pseudo- 
frequencies,"  though  we  should  remember  that  they  don’t  represent  frequencies  of  the  spectrum,  since 
the  spectrum  is  actually  continuous);  these  "pseudo-frequencies"  givp  impulses  when  transformed. 

When  these  impulses  are  convolved  with  the  triangles  due  to  ,  they  give  the  triangle-shaped 

<1> 

sidelobes  on  each  side  of  the  match  point.  Note  that  the  amplitudes  of  the  cosine  terms  in  Eq.  (13)  are 
2,  but  in  the  transform  process  each  cosine  appears  as  2  impulses  whose  amplitudes  are  1,  so  that  the 
amplitudes  of  the  sidelobes  are  1  in  this  specific  code.  Note  also  that  the  sign  of  the  "pseudo- 
frequencies"  determines  the  sign  of  the  sidelobe  (in  this  example,  all  the  sidelobes  are  negative). 

Note  also  that  for  this  example,  the  multiplication  of  S(/)  by  S*(f)  caused  several  e~j7*n  terms 
of  the  spectrum  to  disappear;  here  the  e~J(*  and  terms  of  the  spectrum  disappeared  after 
the  multiplication,  resulting  in  zero  level  sidelobes  at  the  corresponding  locations  of  the  autocorrelation 
function  (see  Fig.  2). 

It  is  clear  that  the  last  term  of  the  spectrum  (generally  e~J{N~  l)J*,  and  here  e~Jl2*)  will  never 
disappear  after  the  multiplication  (since  no  other  term  can  cancel  it),  corresponding  to  the  fact  that  the 
furthest  sidelobe  of  such  code  is  always  +1  or  -1. 

Clearly,  these  observations  will  hold  for  any  phase-coded  signal  with  unity  amplitude  (c.g., 
polyphase  codes  like  Frank  codes),  but  to  any  sidelobe  in  the  autocorrelation,  say  of  g  magnitude,  there 
will  be  a  corresponding  2j?  cos  ( K  •  24  +  9)  “pseudo-frequency"  in  the  power  spectrum.  Generally  g 
can  be  bigger  than  1  but  for  polphodes  g  is  required  to  be  smaller  than  1  (0  is  some  angle  that  depends 
on  the  code). 

To  show  this  process  for  even-length  codes,  examine  the  Barker  code  of  length  4.  It  is  known 
that  there  are  two  possibilities  which  we  designate  as  Barker  Codes  4A  and  40. 

BARKER  CODE  4A 

+  +  +  - 
0  0  0  it' 

The  autocorrelation  function  is  shown  in  Fig.  3: 


i‘(/)~ 

-  w  r 

-l^ll 

1  4-  4-  f~>**  - 

(14a) 

$•(/)- 

■(-£ 

{  4 

J  |l  4-  4-  ~ 

(14b) 

and 

U(/)l3  -  Si/)  $*  (/)  -  |“TJ  j^-j  (4  +  2  cos  2*  -  2  cos  6*j.  (IS) 

Again,  the  autocorrelation  function  corresponds  to  the  "pseudo-frequencies*  of  the  power  spec¬ 
trum  m  Eq.  (IS);  the  match  point  is  4,  the  first  sidelobe  is  +1.  the  second  sidelobe  is  -I,  and  :hc 
cos  40  term  is  missing,  resulting  in  zero  level  at  the  corresponding  point  of  Fig.  3  ( K  -  2). 

Note  that  here,  for  an  even  length  code,  the  signs  of  the  ‘pseudo-frequencies*  cos  2$.  cos  64  are 
opposite,  which  results  in  opposite  sign  sidelobes  in  #(r).  this  property  is  true  for  any  even-length 
(A  -  /J)  Barker  code  (!},  that  means; 
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RESULT  OF 


V  -2  -1H  jo  1  2  V  *(£) 

Fig.  3  —  Autocorrelation  function  of  Barker  code  4A 


R(K)  +  R(N-K)-  0  (16) 

for  any  K  «  1,  2 . tV  -  1,  or  equivalently,  in  |S(/)P  we  will  have  for  any  ±  2  cos  A2$  term  a 

corresponding  =F2  cos  (N  -  K)2<li  term,  such  that  their  signs  are  opposite.  This  also  means  that 
N  jV 

R  K  -  —  -  0  since  the  point  K  -  —  does  not  have  an  ‘image."  Figure  4  shows  the  image  structure 
of  R(r).  The  point  A'  -  1  is  the  "image"  of  A  -  N  -  1,  A'  -  2  is  the  "image"  of  K  -  N  -  2,  etc. 


0  12  3 


R  IK  -  0)  «  N 


"tMAQES" 


N-2  N1  „/! 


ROt  **  — 1»0 


Ftg  4  -  Image  nurture  of  the  autocorrelation  function  of  an  even-length  Barker  tank 


BARKER  CODE  4B 

+  4-  -  + 

0  0  v  0 


6V)-  ~  |^J  [l  +  *“'**  -  e~'*  +  e~^\. 

W>  -  (-^f  **  j~j  (l  +  +  H* 


IS(/)P  —  |-~r,  |  |  U-2  cos  2tfr  +  2  cos  6tfr|. 
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The  correlation  function  is  shown  in  Fig.  5.  Figure  5  is  "similar"  to  Fig.  3.  But  now  in  the  power 
spectrum  (Eq.  (18)),  the  2  cos  20,  2  ccs  60  terms  both  have  changed  signs  when  compared  to  Eq. 
(15),  so  that  Eq.  (16)  is  fulfilled.  This  caused  the  sidelobes  in  Fig.  5  to  change  signs  when  compared  to 
Fig.  3. 


Fi#.  5  -  Autocorrelation  function  of  Barker  code  48 


Note  that  Eq.  (16)  does  not  hold  generally  for  polphodes. 

DEFINITION:  Define  a  O-polphode  as  a  polphode  in  which 

R{K)  +  RHH~  A')-0.  (16a) 

This  is  actually  a  generalization  of  Eq.  (16).  Notice  that  R(K)  can  be  a  complex  number  in  gen* 


As  an  example,  examine  the  generalized  Barker  code  4  (3): 


1  J  -1  J 

0(:  0  re/2  «r  v/2 


5V)- 

i*(/) 


L 

2,V 


|l  +  -  1  •  +  j  • 


and 


-Jl)  e* 

sin  0 

1 

2N\ 

0 

1 

1,  |  7*  1  1 

sin  0 

j 

‘"Wl 

0 

1  <?'*  -  l 


(19a) 

(19b) 


(20) 


corresponding  to  the  values  of  the  autocorrelation  function: 

RiK  -  0)  -  4.  -  1)  -  Rtf  -  2)  -  0.  RIK  -  3)  -  J. 


We  clearly  sec  that  Eq.  (16a)  is  fulfilled,  which  means  that  the  above  code  is  a  G*po!phodc. 
SYNTHESIS  ATTEMPT 

With  the  above  analysis  we  now  try  to  synthesize  the  Barker  type  autocorrelation  function  for 
even-length  (Ar  -  P)  codes.  Barker  and  G- polphode. 
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Barker  Code 


Suppose  there  exists  an  even-iength  Barker  code  xit  x2,  ....  xN.  x,  is  either  +1  or  -1. 
The  autocorrelation  function  of  the  code  must  be  as  illustrated  in  Fig.  6. 


The  match  point  R  (K  -  0)  -  N  appears.  The  nearest  and  furthest  sidelobes  must  appear  with 
magnitude  1  and  opposite  signs  (corresponding  to  cos  20  and  cos  (iV  -  1)  20  terms  of  the  power  spec- 
trum).  In 
AM 
2 


K  -  ~ 


Fig.  6  we  plotted  arbitrarily  one  of  the  two  possibilities  for  these  sidelobes.  At  the  midpoint 
-  0  as  explained  before.  The  dotted  sidelobes  in  the  figure  might  or  might  not  appear. 


But  if  one  dotted  sidclobe  (say  of  index  K)  appears,  there  will  be  a  corresponding  “image“  sidelobe  (of 
index  N  -  K)  with  the  opposite  sign,  as  required  by  Eq.  (16).  Of  course,  there  will  be  another  two 
sidelobes  on  the  other  side  of  the  match  point  (negative  r). 


Now,  from  looking  at  the  desired  R  it)  in  Fig.  6,  we  can  determine  the  structure  of  the  power 
spectrum: 


f* 

sin  0 

4A?} 

0 

(21) 


N  ±  2  cos  20  ±  2  cos  40  ±  , . .  +  0  •  2  cos 


v\ 


:tnust  appear 


±  ...  $  2  cos  (<V  -  2)20  ♦ 2 cos  iff  -  l)20j. 

. 


The  spectrum  Si/)  of  the  code  is  given  *'y  Eq.  (10).  The  magnitude  of  ,Y(/)  must  equal  the 
square  root  of  the  power  spectrum  |A(/)I}  at  every  point  0  (0  was  defined  in  Eq.  (9)  and  represents 
the  frequency  variable).  Specifically,  at  the  A' sampling  points 

0  —  0,  0  -  w/S\  0  «•  2tr/A’,  ....  0  -  fff/fV.  . . . ,  0  »  {S  -  \)#/N 
I/-0.  /-  1  IT.  /-  VT.  ... .  /-  i/r.  ..../-  — ^-1 
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we  must  have: 


|S(/)| 


'iS(/)|2 


IT 

I il  —  i - 

N 


IT 

-  /  — 
/V 


This  is  a  necessary  condition  for  the  existence  of  even-length  Barker  codes,  but  might  not  be  a 
sufficient  condition.  Actually,  Eq.  (22)  gives  us  a  set  of  N  equations  that  must  be  fulfilled. 

Note  that  in  Eq.  (21): 

cos  2  •  0  —  cos  (N  -  1)  2  •  0 
for  &  -  0  cos  4  •  0  -  cos  (N  -  2)  2  ■  0 
cos  6  •  0  —  cos  (N  -  3)  2  •  0 
etc. 

cos  2  ■  rr/N  -  cos  (A/  -  1)  2rr/N 
for  —  w/N  cos  4  rr/N  -  cos  { N  -  2)  2rr/N 
eos  6  t r/S  m  cos  (Af  -  3)  frr/S' 
etc. 


cos  2  •  irr/N  -  cos  (Af  -  I)  2  irr/N 
for  **  irr/N  cos  4  hr/N  •  cos  (N  ~  2)  2  irr/N 
cos  6  iv/N  »  cos  (N  -  3)  2  irt/N 
etc. 

or  generally: 

cos  A  2  irr/N  -  Cos  ( .V  -  A' )  2  ■  hr/N. 

This  means  that  the  power  spectrum  at  the  N  sampling  points  &  -  irr/N  U  -  0, 
(sccEqs.  (16).  (2D); 


iw)i'  |  -  [^j 


*‘2v 


fir 


N  -f  0  4*  0  +  ...■»■  0: 


,v-  1),  is 


i.c..  the  power  spectrum  samples  at  $  «  itt/ N  U  -  0,  1 . N  -  D  must  be  some  constant  N  • 

4, Vs 

times  (sin  d»/*)J  (the  last  term  was  interpreted  ss  the  contribution  oi'  the  basic  code  clement 


length  77  »V). 


Now  the  spectrum  in  those  N  sampling  points  (see  liq.  (10))  is: 


H  "f'wlHMI 


i  v*-  + 


it 

*-*5 


/*»+... +  *  VH 


9 


(25) 
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Denote 


.2  it 


N 


e  "  -  W 

(this  is  the  known  basic  phasor  of  DFT  where  WN  =  1). 


(26) 


S(f) 


T 

sin  *li 

.1 

in 

2  N 

e 

* 

l 

in 

•  [l  ■  +  W'eJ*2  +  W2ieJ*3  +  . . .  +  Frw“I>'  eJ*N\  (27) 


il> - 

N 


i/i - 

N 


and  requiring  (22)  results  in  )V  equations: 


/-  0: 


ll  •  eHl  +  1  •  eJ*2  +  1  •  eJ*'+  ...+  1  •  JN  -  / 

1  .  J*\  +  w  .  eJ*i  +  w1  eJh  +  ...  +  WN~X  eJ*N\ -y/N  -  / 
r  •  eJ*x  +  W2  e1*2  +  W*  eHi  +  ...  +  eJ*N\ - VF  -  / 


(28.1) 

(28.2) 

(28.3) 


W-  1: 


1  •  e7*1  +  WN~X  e1*1  +  WN~2  eJ*}  +  . . .  +  WeJ*N  -  y/N  —  / 


and  in  matrix  notation: 


ii  i 

1  w  w2  . 

1  w2  w* 

1 

WN- 1 
WN-  2 

eJ*> 

eJ*> 

'ten 

l/<*2 

!/<*) 

1  WN~2  . 

W 

eJ*s 

I/<*N 

(28.N) 


(29) 


Going  from  Eq.  (28)  to  Eq.  (29),  we  had  to  take  care  of  the  absolute  value  in  the  left  side  of 
(28),  by  placing  some  'in known  phases  a,  in  the  right  sido  of  (29)  for  each  element  whose  magnitude 
should  be  exactly  /  -  \N. 

We  can  write  Eq.  (29)  as: 

A  X  -  V  (29a) 

where  A  is  the  known  DFT  matrix  (N  x  /V  matrix),  which  is  nonsingular  with  dot  A  &  0. 
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The  phasor  is  on  the  unit  circle  (see  Fig.  7). 


*1 

e^\  ' 

*2 

eJ*: 

*3 

>4 

c  • 

% 

is  our  unknown  vector,  which  represents  the  required  Barker  code  (x,  -  ±  1,  ot  is  either  0  or  ir). 


‘lSL 

l/a2 

V2 

//«3 

- 

v3 

l/<*N 

Vs 

is  a  vector  whose 


elements  have  magnitude  /  -  \fN,  with  unknown  phases  a,. 


COMMENTS 

(1)  Equations  (29)  are  exact  necessary  conditions. 

(2)  Equations  (29)  hold  only  Ter  even-length  codes  N  -  /J.  A  similar  analysis  for  Barker  codes  5,  7, 
11,  and  13  shows  that  the  spectrum  samples  (of  the  sequence)  are  not  required  to  have  a  constant 
magnitude. 


II 
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G-polphode 

A  similar  analysis  for  the  G-polphode,  where  Eq.  (16a)  is  fulfilled,  will  give  us  the  same  result. 
For  any  time  sidelobe  R  ( K )  there  will  be  a  corresponding  R{N  -  K)  time  sidelobe  such  that, 

R(K)  +R*(N-K)~  0 

or  equivalently:  Real  (/?(£)]  +  Real  [A  (A  -  A)]  «=  0 

IJR (K)]  -  IJR(N-K))-  0. 

This  means  that  the  pair  of  sidelobes  R  (A')  and  R(N  —  K)  contribute  to  the  power  spectrum: 

2  Real  [R  (AT)}  (cos  K  ■  2«/r  -  cos  (N  -  K) 2^}  +  2 Im[R  (A')l  (sin  K  •  20  +  sin  (N  -  K)2p}. 


This  contribution  of  the  pair  goes  to  zero  for  the  N sampling  points  p  **  in/  N,  since 

cos  K  ■  2—r  -  cos  (N-  A)2-~ 

N  N 

sin  K  ‘2—  **  -  sin(yV—  K)2  ~ 

N  N 


thus  resulting  in 


I  S(/)l2 

<l> 


hr 

N 


T1  s[nJ/_ 
4/^  1 1> 

* 


•  [N  +  0  +  0  +  0  +  ...  +  0) 

In 


as  before. 


So  Eqs.  (29)  and  (29a)  hold  also  for  G-polphodes,  but  the  code  elements  can  be  any  complex 
number  with  unity  magnitude  U,l  «  I. 

Thus,  from  now  on  we  can  proceed  with  a  a  sequence  qj  numbers  X,  (real  for  Barker  and  complex 
for  G-polphode)  which  when  DFT  transformed  (Eq.  (29)).  gives  a  vector  with  constant  magnitude  ele¬ 
ments  /  -  y/W. 


We  will  examine  first  Barker  codes, 


BARKER  CODE  STRUCTURE  AND  PROPERTIES 

To  derive  several  properties  of  an  even-length  Barker  code  (if  il  exists),  we  write  the  mapping  Eq. 
(29)  m  a  convenient  form: 


A|  + 

1 

•Vj  +  l  •  .Vj 

+ 1 

•  a4  + 

...  +  l 

•  a -ft  -  t[n± 

01.1) 

At  + 

U> 

•  A,  +  IP*’  • 

•VJ 

+  w*  • 

A«  +  . . 

.+  IP*-'  •  jr*  - 

'  ‘bi 

01.2) 

•V|  + 

•  Xj  +  W4 

‘  Aj 

+  w* 

•  A<  +  . , 

. .  +  tP'V~3  •  XN  -  l{ Oi 

01.3) 
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1  •  x{~  1  ■  Xj,  4  1  •  X3  -  1  •  A4  4  . . .  -  1  •  XN  -  l[aN/ j+  , 


(31.N/2  +  1) 


1  A,  4  WN~2  •  X2  +  WN~*  •  *3  4  WN -6  •  A4  +  . . .+  W1  •  Xh  -  //ttjy-1 
1  •  A,  +  ■  Aj  4  WN~2  •  Aj  4  H"*"3  •  A4  4  ...  4  W  •  A*-  //«w 


(31.N-1) 

(31.N) 


from  £17.  (3/..0:  since  a,  is  real  (±  1),  oel  must  be  0  or  tr,  so  that: 

At  4  A2  +  A3  +  . . .  4  XN  -  ±  /, 
i.e., 

number  of  pluses  -  number  of  minuses  -  ±  /.  (32) 

But  since  their  sum  is  N  -  /J;  then: 

1 2  +  /  i2-l 

CASE  1:  if  number  of  pluses  -  — - —  then  number  of  minuses  -  — y- 

(e.p  ,  Barkers  4  4-4  and  4  +  4  -) 

I2- 1  1 2  4  / 

CASE  2:  if  number  of  pluses  -  -y—  then  number  of  minuses  -  — - — 

(e.g.,  Barkers  4 - and  -  + - ) 

For  simplicity  we’tt  discuss  only  Case  1  in  the  following  few  paragraphs  (Case  2  is  the  "opposite"  case). 

Note  that  the  difference  be.wccn  the  number  of  pluses  and  minuses  gets  larger  as  the  code  length 
increases,  which  is  not  toe  case  in  PN  b!".ary  sequences. 

From  Eq.  (3LN/2+I):  agrln  a^,  must  bo  ^  or  ir,  and  : 

A|  -  A*  4  X\  -  A«  4  . . .  4  A/v_|  -  Xtf  —  ±  /.  (33) 

Odd  pluses  and  oven  minuses  contribute  positive  numbers  in  Eq.  (33),  while  even  pluses  and  odd 
minuses  contribute  negative  numbers. 


Denote: 


m  «  number  of  odd  pluses,  then— - - m  «  number  of  even  pluses 

I2  -  / 

»  -  number  of  odd  minuses,  then— y-  -  «  -  number  of  even  minuses 


From  Eq.  (33): 


f  /3  4  / 

|  - - -  m  4  n 


2  (m  -  n)  -  /  "•  ±  / 
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We  have  two  possibilities: 


(A)  m  =  n,  but  since  m  +  n  =  —  (number  of  all  the  odd  element.^,  we  get: 

/2 

number  of  odd  pluses  -  number  of  odd  minuses  -  — 
(e.g.,  Barker  +  +  -  +  ) 


(B)  m  -  n  —  I,  this  implies  similarly  that: 

[2 

number  of  even  pluses  -  number  of  even  minuses  -  —  (34.2) 

(e.g.,  Barker  +  +  +  -  ) 

From  Eqs.  (31.2)  and  (31. N):  each  weight  of  the  real  code  elements  (x,  -  ±1)  in  Eq.  (31.2)  is  the 
complex  conjugate  of  the  corresponding  weight  in  (31. N),  e.g.,  W*  -  WN'\  (IV2)*  -  WN~2,  etc.,  so 
that  Ifct}  must  be  the  complex  conjugate  of  l/aN,  or: 


<*n  “  -«2- 


(35.1) 


Similarly: 


«yv-l  -  -a3 

«W-J  “  ~«4 


(35.2) 

(35.3) 


aN/)  +  2  -  -a#/?- 


(35.4) 


These  equations  say  that  for  real  codes,  Eqs.  (29)  take  the  form: 


//0  or  ir 


i 

*”  >tSi 


NRL  REPORT  8586 


From  Eqs.  (31.1)  and  (31.N/2+1):  by  adding  and  subtracting,  we  get  (taking  into  account  Eq. 
(36)): 


0 

2(jcj  +  x3  +  x$  +  x7  +  . . .+  xyv_i)  -  or  (37a) 

/ 

/ 

2(x2  +  x4  +  x6  +  x8  +  . . ,  +  xN)  -  or  (37b) 

0 


which  are  another  form  of  Eq.  (34). 

Further  properties  of  Barker  codes  can  be  derived  if  one  can  follow  the  requirements  logically.  As 
an  example,  consider  the  UV/4+1)  which  is  a  member  of  Eq.  31,  and  its  conjugate.  With  property  (36) 
in  mind  (note  that  these  two  equations  give  ±  90°  shift  in  the  weight  from  each  code  element  to 
another),  we  get: 

1  •  X\-p Cj  -  +  jt*4  +  1  •  X$-JX(,-  1  •  JC7  +  Jxt  +  . ..  -  l/fl  ,  (38.1) 

and 

1  •  xt  +  jxj  -  1  •  Xj- JX4  +  1  •  Xj  +  jSx$- 1  •  x2  -  jSX|  +  . . .  -  I/-8 .  (38.2) 

By  adding  and  subtracting  we  get 

2(xl-X}  +  XJ-X7  +  X9~Xlt  +  Xij-Xjj  +  ...)  -  //g+  //-] 8  .  (39.1) 

and 

2y(-xj  +  X4~X4  +  X|-X|o4-Xjj-Xi4  4-Xt6-...)-  //§-  //-) 8  ,  (39.2) 

or 

(Xj-Xj  +  XJ-X7 +...)-/ cos^  .  (40.1) 

and 

(-Xj  +  X4-X4  +  xt- ...)-/  sin  /3.  (40.2) 

Equations  (40)  can  be  fulfilled  simultanously  for  a  few  possibilities  of  the  angle  fi,  since  their  left 
side  is  on  integer  (with  plus  or  minus  sign).  Actually,  if  It*  Sp(not  multiple  of  S),  the  only  values  for 
8  are  0,  ±  90®,  ±  180°,  which  result  in  an  integer  on  the  right  side  of  Eqs.  (40.1)  and  (40.2).  If 
I  -  5 p  (multiple  of  S),  there  are  other  possibilities  to  get  an  integer  in  the  right  side,  since 
cos  8  -  3/5  or  cos  8  -  4/5  results  in  sin  8  "  4/5  or  sin  8  "*  3/5,  which  means  we  have  another  ■fam¬ 
ily'’  of  possibilities  that  can  fulfill  Eqs.  (40).  Actually,  they  arc  all  the  possible  combinations  of  ±  3/5, 
±  4/5  for  the  cos  8>  sin  8< of  Eqs.  (40). 

Another  important  observation  is  derived  by  adding  all  the  equations  of  (31).  Then  in  the  left 
side,  all  the  code  elements,  except  .V|,  will  cancel  (because  the  weights  are  uniformly  distributed  pha- 
sors  in  the  unity  circle  of  the  complex  plane),  resulting  in: 

I*  X)  ■■  I/a  |  +  //«2  +  //at)  +  . . .  +  la #  , 
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and  since  we  can  assume  x{  —  1, 

I2  »  lj_ 0  +  //ft;  +  1/ft;  -f  . . .  +  //0  or  7T  +  . . .  +  //— «3  +  //"ft;  .  (41) 

which  means  that  all  the  l2  phasors  in  the  right  side  of  (31)  or  (36),  whose  magnitudes  are  4  and  which 
appear  in  pairs  of  complex  conjugates,  must  sum  to  I2.  This  means  also  that  one  possible  choice  of  the 
phasor’s  vector  in  the  right  side  of  (31)  is  the  code  itself  X  times  L  In  such  a  case,  the  right  side  of 
(31)  is: 

Ixi  +  bci  +  . . .  +  bcN  -  /Of)  +  x;  +  . . .  +  xN)  —  I  •  l  -  I2  (42) 


i 

i 

i 


as  required  by  (41). 

All  the  above  properties  ((31)  through  (42))  can  be  utilized  to  reduce  the  search  for  even-length 
Barker  codes. 

PHYSICAL  INTERPRETATION  FOR  BARKER  CODES 

We  can  examine  now  the  physical  meaning  of  Eq.  (36),  as  illustrated  in  Fig.  8.  We  need  to  input 
the  real  code  x,(il)  to  a  DFT  system,  such  that  we  get  a  constant  amplitude  /  in  the  output,  while  the 
phases  of  the  output  must  fulfill  some  constraints. 


TIME  DOMAIN  FREQUENCY  DOMAIN 

FI*.  I  —  Ptiytkii  meaning  of  Eq.  (16):  the  DFT  of  the  real  sequence  X,  fives  constant  munitude  phasors 
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Now  we’ll  see  how  the  Barker  codes  for  N  -  4  (/  «•  2)  are  derived  by  the  above  analysis  (see  Fig. 


fy  _  ^-7  2*/V  _  e~J  2ir/4-  -> 


1111 

1111 

1  w  W1  wl 

1  -7-1  J 

1  W1  w*  w6 

1-1  1-1 

1  Wl  ff*  H* 

1  7-1-/ 

We  need: 


*1 

2/0 

*2 

2/a 

2/ir  or  0 

x* 

l2 Lzi 

We  see  from  Fig.  9  that  C\  only  can  create  the  dc-frequency  term  2/0,  C2  only  can  create  the 
fundamental  frequency  2 /a,  etc.  Thus  the  required  code  Cis  a  linear  comotnation  of  Cj,  Cj,  Cj,  C4 
( in  the  time  domain).  Ifwe  can  find  a  code  C  all  of  whose  elements  are  of  unity  magnitude,  then  it  is 
the  required  code  (note  that  C3  has  two  possibilities). 


C,: 

1/2 

1/2 

1/2 

1/2 

C2: 

1/2 /a 

-/•  1/2/a 

-1-1/2/a 

+J  •  1/2  /a 

Cj:  (g) 

1/2 

-1/2 

1/2 

-1/2 

-1/2 

1/2 

-1/2 

1/2 

C4: 

1/2 /^a 

-11/2/^ 

+/•  1/2 /-a 

C. 

*2 

Xi 

*4 

We  have  only  one  parameter  (a)  to  choose  in  order  to  have  the  required  code,  all  of  whose  elements 
must  have  unity  magnitude.  Wc  sec  that  if  Cj®  is  examined,  a  must  be  +90"  or  -90®  (from  the  first 
column,  in  order  to  have  xj  ~  1),  so  the  code  is: 


C:  xt  -*  1,  x}  -  1.  xj  »  1,  x4  -  -1  for  a  «*  90® 
C:  x(  «  1,  Xj  *»  -1,  xj  «•  1.  x«  «■  1  for  a  -  -90® 
and  if  Cj(8)  is  examined,  a  must  be  0®  or  180®: 

C:  x,  -  1.  x,-  1.  x,-  -1.  x4~  I  fora-  0“ 
C:  x,  »-!.  Xj~  1.  x,-  1.  x«-  1  fora  -  180® 


I"  All  the  above  codes  C  are  legitimate  Baker  codes  which  fulfill  all  the  requirements.  Here  for  S  m 

&  4,  wc  had  only  parameter  a  to  choose,  but  when  A/  is  large,  wc  have  many  parameters  to  choose,  such 

that  all  the  elements  in  C  will  add  up  to  unity. 


A  pictorial  interpretation  of  the  requirement  established  by  (36)  is  illustrated  in  Fig.  10  (only  for 
Barker  codes)  for  three  elements  of  the  vector  matrix  described  by  that  equation. 


FI*.  10  —  Pictorial  Imwpteuiioo  of  Eq.  06) 


LINEAR  ALGEBRA  POINT  OF  VIEW 

We  now  analyze  our  problem  for  either  Barker  codes  or  G-polphodes.  Equation  (29),  which  is  a 
necessary  condition  for  both  of  them,  can  be  written  as: 
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Define: 


b-7a 


which  is  a  modified  DFT  matrix  (each  element  in  A  is  divided  by  /  in  order  to  get  B. 
Then: 


B  X  -  Y, 


where  Y  is  the  vector: 


Equation  (44)  requires  the  code  vector  X  to  map  to  vector  Y  (unity  magnitude  elements)  through 
the  modified  DFT  matrix  B. 

This  can  happen  in  two  ways: 

1.  The  vector  Y  is  some  scalar  A  (might  be  complex)  times  X.  Then: 

BX  -  XI  (45) 

We  will  call  this  case  an  eigenvector  mapping  code  (we  have  mentioned  this  possibility  for  Barker 
codes  after  (41)) 

2.  Y  *  A  g  (46) 

We  will  call  this  case  a  noncigcnvcotor  mapping  code. 

In  order  to  investigate  the  eigenvector  mapping  case,  we  will  use  some  properties  of  the  matrix  B 
(over  the  complex  field). 

Writing  (44)  in  detail,  we  get: 


1 

1 

1 

1 

7 

/ 

w 

Wi 

1 

1 

t 

wl 

w « 

WN-1 

1 

/ 

1 

ws~* 

WJ 

/ 

/ 

/ 

M*-' 

w 

/ 

/ 

1 

X\ 

Xj 

r> 

Y*. , 

r» 

where  U,|  -  1,  1  T,|  -  1  /  *»  1,  2 . Af. 
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PROPERTIES  OF B 

a.  The  columns  (or  rows)  of  B  are  orthonormal: 

(±i)r (JS*)-  i  \tJj  (48) 

where  Vt  denotes  the  /,h  column  vector  (note  that  this  is  the  definition  over  the  complex  field,  as  a  gen¬ 
eralization  of  orthonormality  over  the  real  field). 

b.  B  is  symmetric: 

B  -  BT.  (49) 

Also,  its  rows  (except  the  first  and  (jV/2  +  l),h)  are  pairs  of  complex  conjugates,  e.g.,  the  /Vth  row 
is  the  complex  conjugate  of  the  2“*  row,  the  (N  -  1)*  row  is  the  complex  conjugate  of  the  3rd  row,  etc. 


c.  B  is  a  unitary  matrix  (this  is  the  complex  generalization  of  an  orthogonal  matrix  over  the  real 
field,  where  A  is  an  orthogonal  matrix  if  AA  r-  /),  which  is  defined  by: 


B(B*)T -  /. 

($0) 

or  equivalently: 

( B*)t . 

(50a) 

and  in  our  case,  due  to  (49): 

-  8 * 

(51) 

From  (47),  (51): 

I-  **  K 

(52) 

d.  Idet  B 1  -  l 

(53) 

for  any  unitary  matrix  (see  (4),  p.  112),  which  means  that  B  is  a  nonsingular  matrix  of  rank  N. 

e.  All  the  A' eigenvalues  of  B  (as  a  unitary  matrix)  are  of  unity  magnitude  (see  (4)  p.  15$,  prob. 

22); 

lx,l-l  <-1,2 . S'.  ($4) 


It  can  be  verified  that  in  our  case,  at  least  Xt  •  1,  Xj  -  -I  arc  i  gen  values  of  B,  possibly  with 
some  multiplicity.  To  show  this: 
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Adding  all  the  rows  of  (55)  to  the  last  row  we  get: 


Ift-A/I- 


|I_X  1  ...  i 
1  /  /  / 


/  —  A  —A  ...  —A 


(56) 


This  last  determinant  is  zero  for  A!  —  1  and  A  2  ■*  -1  (since  for  both  of  them  we  get  two  propor¬ 
tional  rows  in  the  determinant). 

f.  ft  as  a  unitary  matrix,  maps  any  vector  X  to  vector  Y,  such  that  their  energies  arc  the  same 
(mathematicians  call  this  property  preservation  of  length)  i.e.: 


*1  X\  +  Xi  xl  +  ...  +  XN  Xs  -  Yx  y;  + 
Yi  y:*  +  ...+  y;. 


(57) 


ft- 


Note,  however,  that  if  |x,|  -  1  (unity  magnitude  code)  Y,  generally  are  not  necessarily  of  unity 
magnitude.  Our  problem  is  to  find  that  lx,  I  ■»  1  that  will  map  to  i  y(  I  -  l ,  and,  of  course,  it  is  possible 
from  an  energy  point  of  view. 

As  an  example,  check  the  case  /  —  2  (N  »  4): 

1  i  1  1 

2  2  2  2 

I  ^  d  i 
2222 

1  1  1  1 

2  “2  2  ~2 

i  id  d 
2222 

If  the  code  is  an  eigenvector  of  ft  then,  BX  “  A  A' 

The  eigenvalues  are,  Aj  *  +1.  Aj  -  + !.  A>«  -1,  A4  -  -j 
so  that: 

del  (ft  -  A/1  -  (A  -  1)  U  -  1)  (A  +  I)  (A  +  j).  (58) 

Note  that  indeed  U,  i  »  1  and  idet  ft  I  -  i/i  *  1  (when  substituting  A  *•  0  in  (58)). 

The  eigenvectors  are: 

l.  for  A,  -  Aj  —  1  we  have  two  eigenvectors: 


(ft  -  /)  V  -  0  ~ 


0  0  0  0! 

0  0  0  0 

10-1-2 
0  1  0-1 
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1 

1 

1 

0 

Jfl- 

-1 

and  Vi  - 

1 

1 

0 

'V': 


■  V-1' 

ft  v 


V\  is  a  Barker  code,  while  V2  is  not.  Actually,  any  linear  combination  of  V\  and  V2  that  has 
constant  amplitude  is  also  a  good  solution  (in  our  case  only  V\  and  -V\  are  Barker  codes). 


2.  for  Xj  -  -1, 


(B  +  /)  V  -  0 


0  0  0 
2  1  0 
1  0  1 


1  j  -l  -j  +  2 


0. 


Kj  is  the  eigenvector  Barker  code  (of  course,  -V)  is  also  a  good  solution). 


3.  for  X«  -  ~j. 


(B  +  Jl)  jr-o~ 


f4- 


0  111 
0  10  1 
l  0  0  C 
0  0  0  0 

0 

i 

0  • 

-1 


V  -0. 


The  eigenvector  V4  is  not  a  Barker  code.  Now  we  will  prove  that  for  the  eigenvector  mapping 
BX  -  kX,  only  X  -  +1,  X  -  -1  can  give  us  a  legitimate  Barker  code  or  G-polphede  (where  U‘(|  —  1). 


For  an  eigenvector  ntapping  we  require 

[i  I  i  ...I 


III  1 

l  W  jf* 

,V| 

*1 

Ai 

i  ir5  w* 

A) 

Aj 

7  /  i  ... 

. 

-  X 

- 

• 

% 

Ajf 

1  ws-'  w*-}  ... 
7  I  t 


(59) 
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From  the  first  row: 

-y  (*!  +  A 2  +  .  . .  +  Xjy)  =  XJCj.  (60) 

From  summing  all  the  Eqs.  in  (59)  we  get: 

(-y  X)  +  -y  X|  +  . .  .  +  -y  .\)  +  0  •  X2  +  0  •  *3  +  •  •  •  +  0  •  X/y  =  X(xt  +  X2  +  .  .  .  +  XN ).  (61) 

From  (60)  and  (61): 

I2  ■  X|  =  X  •  X  lx\,  (62) 

t 

or 

X(  -  X2  Xi-  (62a) 

Equation  (62a)  can  be  fulfilled  only  if: 

a.  X!  =  0,  which  will  not  give  a  Barker  code  or  G-polphode  (requires  |x,  |  -  1), 

b.  X2 »  1  or, 

X  -  ±  1.  (63) 

which  might  give  a  Barker  code  or  a  G-polphode. 

Thus,  an  eigenvector  code  can  be  achieved  only  for  the  eigenvalues  X  -  ±  1. 

The  other  complex  eigenvalues  jx,i  -  1  will  not  give  a  desired  code  (we  saw  it  in  the  example  for 
/  »  2,  where  X4  ~  -j  did  not  j:*e  a  Barker  code,  and,  indeed,  the  first  element  of  the  eigenvector  Vx 
was  Xi  =  0). 

Now  we  prove  that  an  eigenvector  mapping  does  not  have  a  solution  for  a  Barker  code 
(x,  »  ±  1 )  for  /  > 


If  A'  -  (xj,  x2 . x/y)ris  r-a/,  then  the  eigenvector  possibilities  arc: 

for  X  ™  v  BX  -  \  •  X  and  (64a) 

for  X  —  1 :  BX  -  - 1  •  X,  (64b) 


Before  proceeding  with  the  proof,  it  will  help  to  observe  the  ease  /  *  2. 


a.  BX  -  1  X  gives  the  eigenvector  Barker  codes 

F,-(++  +)r 

r 


b.  BX  -  -l  •  Af  gives  the  eigenvector  Barker  codes 

V2  “  (—  +  4-  +)r 

—if,  —  (+ - )r- 

But  note  that  «■  (+  ~  +  +) r  and  — Jfj  -  (-  +  — )r  arc  not  eigenvectors  of  B ,  though  they 
are  Barker  codes,  which  are  obviously  "symmetrical"  to  the  above  V\  and  -V\.  For  example. 
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1  i  1  1 

2  2  2  2 


1  zi  _i  1 

'  l1 

1 

2  2  2  2 

-1 

j 

1111 

* 

1 

1 

2  2  2  2 

1 

~j 

1  1  _I  zi 

2  2  2  2 


is  not  an  eigenvector  though  it  is  a  Barker  code.  Similarly,  _K4-(+  +  +  -)r  and 
-_P4  =  ( - +) 7  (which  are  symmetrical  to  the  above  _K2  and  -Jfj)  are  Barker  codes  but  not  eigen¬ 

vectors  of  B.  This  happens  because  of  ihe  general  requirement  that  a  real  eigenvector  must  obey  the 
follc'ving  structure  (for  X  «  ±1): 


*\ 

*1 

*2 

Xj 

*} 

-  X 

■ 

■XN~\ 

X*3 

*N 

x*2. 

(65) 


This  was  explained  in  Eq.  (36)  and  Fig.  10,  for  a  real  code  X.  But  if  x,  is  real  (±1),  then  x,  -  x’, 
so  that  Eq.  (651  requires: 

%  "  xj  -  Xj, 

*N~l  "  X)  "  X), 

xN_j  -  xj  -  x4,  etc.  (66) 


We  see  thut  JK,.  ~V ,,  V},  -Vj  above  fulfill  this  requirement  (.vj  -  x4),  and  th:  rcfore  can  be  real 
eigenvectors.  On  the  other  hand.  jfj.  -Vy  Jf4,  -Jf4  do  not  fulfill  (66)  (since  x4*  ~x2),  and  there¬ 
fore  cannot  be  eigenvectors. 

Now  to  proceed  with  the  proof,  the  next  candidate  for  our  problem  is  /  *•  4  (N  -  16). 

According  to  the  above  analysts,  for  the  real  eigenvector  mapping,  the  code  structure  must  fulfill 
Eq,  (66).  Thus,  the  eigenvector  code  must  be: 

CODE:  X|  x?  xj  x4  x5  xf  ...  x6  xj  x4  .v,  xj 

ELEMENT  NO,  1  2  3  4  5  6  ...  N-4  N-3  N-2  N-l  N 

whore  x,  is  cither  1  or  -  I  (note  that  the  above  £],.  -P’1.  £2.  -£2,  fulfill  this  structure). 

To  show  that  this  is  impossible  for  /  >  2,  wc  return  to  the  time  domain  autocorrelation  process 
by  steps. 

RRST  STEP 

.  .  .  X»  X4  Xj  Xj 

X,  .Vj  Xj  . . . 

xj.  Xj  can  be  A 1,  so  that  R{N  -  1)  «•  ±1. 
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SECOND  STEP 


*5 

A'4 

*3 

*2 

*1 

*2  *3  *4  *5  ■  •  • 

R  (N  —  2)  =  X2  ■  x2  +  *1  ■  *3  “  1  +  *1  •  *3- 

Since  R(K)  is  allowed  to  be  0  or  ±  1,  it  follows  that  x2  —  -xt 

THIRD  STEP 

•••  *6 

*5 

*4 

-*1  x2 

*1 

X2  —X!  X4  X$  ... 

R  (N  -  3)  -  X[  (x4  —  2x2)  so  that  x4  ■»  x2. 

FORTH  STEP 

•••  *6 

*5 

*2 

-X,  x2 

At 

*2 

R  (N 

-xt  x2  xs  .  .  . 

-  4)  -  x2  ■  x2  +  X[  •  X]  +  x2  •  x2  +  X|  x5  -  3  +  x( 

No  X|,  x5  (which  are  ±  1)  can  give  the  desired  autocorrelation  function  (0  or  ±  1),  thus  proving  that  no 
real  eigenvector  code  exists  for  /  >  2. 


By  now,  wc  see  (hat  the  remaining  possibilities  to  meet: 
BX-Y 

lx,!-l.  I K, I  -  1 


are: 


1.  Barker  code  (real),  x,  «  ±1: 


B 


*t 

y. 

*2 

y2 

AJ 

m 

Yi 

%-l 

y; 

XN 

Yl 

(68) 


(68.1) 


where  Y  &  k  X  (not  an  eigenvector)  for  /  >  2.  This  possibility  hus  to  meet  properties  (32)  through 
(40). 

2.  G-polphode: 

a .  8  X~  ±1  -  AT eigenvector  mupping.  (Sec  the  appendix  for  further  properties  in  this  case.)  (68.2) 

b.  noncigcn vector  mapping,  Y  &  k,  X. 


As  an  example  of  possibility  2.b,  consider  the  specific  polphodes  that  arc  given  by  the  generalized 
Barker  codes  13j,  These  arc  derived  from  a  ’father"  Barker  code  Xs  by: 

.v,  -  <*,)*  e'1-'*.  (69) 

where  6  is  some  angle  2w/ P  (Pis  an  integer). 
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It  is  actually  the  addition  to  each  element  (x,)  g,  of  a  progressing  phase  step  (6  can  be  further  gen¬ 
eralized).  This  modification  does  not  change  the  envelope  of  the  autocorrelation  function. 

Thus,  all  the  codes  defined  by  (69)  form  polphodes,  some  of  them  are  G-polphodes. 

Examples  for  JV  =  4: 


1 

1 

1 

0-90’ 

j 

Kb - 

1 

— *•  T- 

-1 

-1 

J 

Bi¬ 


ll  ±  ± 

2  2  2  2 

1  zL  zl  1 

2  2  2  2 

1  l  l_l 

2  2  2  2 

2  2  2  2 

We  sec  that  X  is  not  an  eigenvector,  but  it  is  a  G-polphode  since;  R  (3)  ->  j,  R  (1)  -  J,  as 
required  by  (16a). 

2. 


1 

J 

j 

1 

-1 

-7 

J 

i 

1 

1 

1 

*-90* 

7 

la- 

-1 

-*•  X  - 

1 

1 

-7 

1 

1 

7 

1 

1 

1 

-7 

-1 

BX-n 


X  is  not  an  eigenvector,  but  it  is  a  G-polphode  since,  R  (3)  -  J,  R  (1)  -  -J. 


Xu  - 


1 

0-45* 

1  /4S_° 

1  /90; 

1  /— 45 0 

BX  - 


1 

1 

1 

-1 


1  I 

2  2 

i.  zl 
2  2 

1  1 

2  ~  2 

1  1 

2  2 


1 

2 


1 

2 

zi 

2 

1  I 

2  "2 

_1  zl 
2  2 


1  /45° 

1/90; 

IMS 


y  +  y  j  +  cos  45 
y  ~y;  +  sin45°* 
y  +  y  J  -  cos  45 

2  “I^sln4S 


(70) 
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X  is  a  polphode,  but  not  a  G-polphode  since  the  autocorrelation  funtion  is. 
Riff)  =  4,  i?(l)  =>  I  /45°,  Ri 2)  =  0,  /? (3)  =  1  /— 45 °,  and  Ril)  +  R*( 3)  ^  0,  in  contrast  to  Eq. 
(16a).  Note  also  that  the  right  side  of  Eq.  (70)  does  not  have  constant  amplitude  elements.  This  last 
example  shows  that  there  might  be  polphodes  that  are  not  G-polphodes,  thus  our  analysis  does  not 
cover  them. 

At  this  point,  we  review  our  results  as  shown  in  Fig.  11.  A  question  mark  denotes  codes  that 
were  not  investigated  in  this  paper. 


/ 


COOES 


\ 


?  BARKER  TYPE  AUTOCORRELATION  FUNCTION 


REAL  (BARKER) 

Xj  ■  1 1 

/  \ 


OOP  LENGTH 

BARKERS 
3.5. 


EVEN  LENGTH 


KERS  /  \ 

7.  11,13  /  \ 


N“2  N-F  NP./2 
[♦-]  f  \  IMPOSSIBLE 

(-1  /  \ 

N«4  J>  2 

■/*2 

[♦♦-♦]  NON-EIGENVECTOR, 
r  -I  PROPERTIES 
+  (J2)  THROUGH  (401 


\ 


POlPHOOES 

K!" 

/  \ 

OOP  LENGTH  EVEN  LENGTH 

~r'/  \ 

N#l2  N-i2 

“  /~\ 


O-POLPHOOES  NOT  G-POLPHOOES 


/  V 

BX«tl*X  BX-YPXX 

EIGENVECTOR  CODE  NON -EIGENVECTOR 


Fig.  1 1  -  Roview  of  results;  a  question  mark  denotes  codes  (hat  were  not  investigated  in  the  paper 


Z-TRANSFORM  INTERPRETATION 

Further  insight  into  the  problem  of  generating  a  code  is  achieved  by  using  the  Z-transform.  Basi¬ 
cally,  we  need  a  sequence  x,(U,l  -  1,  i  -  1,  2,  ... .  N,  where  A'  -  /2)  such  that  its  DFT  will  have 
constant  magnitude. 

The  DFT  of  a  sequence  is  given  by  N  sampling  points  of  the  Z-transform.  The  sampling  points 
are  uniformly  distributed  on  the  unity  circle  of  the  Z-planq  (see  Fig.  12). 
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SAMPLING  POINTS. 


Fig.  12  -  The  DFT  of  a  sequence  is  given  by  N  sampling  points  (uniformly 
distributed  on  the  unity  circle)  of  the  Z-transform 


The  Z-transform  of  the  sequence  x\,  xj,  ,  xN  is: 

X(Z)  -  x»  +  x2Z-‘  +  xjZ"2  +  . . .  +  xnZ-(n-u 

x,  Z*-'  +  x2  ZN~ 2  4-  x,  Z*"3  +  . . .  +  x/v_,  Z  +  xN 
zn-\ 


Then: 


[DFT  of  x,]  -  XU) 


-  X(.K) 


z  -  WK 


(71) 


(72) 


where  W  -  e~^N,  K  -  0.  1 . N  -  1. 

We  see  in  (71),  that  XU)  has  N  -  1  poles  at  the  origin  (Z  -  0),  and  N  -  1  zeroes  that  depend 
on  the  sequence  x,. 

If  the  x/s  are  real  (±  1),  the  roots  of  the  polynomial  in  (71)  are  either  real  or  complex  conjugates 
in  pairs. 


Since  N  is  even,  N  -  1  is  odd,  so  that  out  of  the  N  -  1  zeroes  of  X(Z)  there  will  be  an  even 
number  of  complex  conjugate  zeroes  and  an  odd  number  of  real  zeroes. 


MV' 


Thus  X(Z)  for  a  real  sequence  x,  can  be  factored  to  the  form: 

X(Z)  -  (Z  -  Z,)(Z  -  Zj)(Z  -  Z3)  . . .  (Z  -  Z4)(Z  -  Z\  . . .  (Z  -  Z5)  (Z  -  Z\)  (73) 


odd  number  of  real  zeroes  pairs  of  complex  conjugate  zeroes 

Since  wc  are  interested  in  the  magnitude  of  the  DFT  of  the  sequence  at  Z  -  WK  -  e~jiwK,N 
where  K  -  0,  ....  N  -  1,  we  can  ignore  the  (N  -  1)  poles  at  the  origin  (they  do  not  affect  the  magni¬ 
tude  of  X(K)). 
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As  an  example,  examine  a  Barker  code  of  length  4: 

*1—1  X2“l  1  *4  -  1 

-3  Z3  +  Z2  -  Z  +  1  (Z  -  Z/)(Z  -  Z2)(Z  —  Z3) 


X(Z)  -  1  +  Z~l-  Z~2  +  Z- 

Carrying  out  the  factorization  we  get: 

Z,  =  -1.84,  Z2  *  0.42  +  y'0.6,  Z3  =  0.42  -  y0.6  -  Z\. 


(74) 


Those  values  are  calculated  approximately  for  the  sake  of  illustration  (see  Fig.  13). 


Fig.  13  -  Poles  and  zeroes  of  the  2-transform  for  Barker  code  4, 

A,  B,  C,  D  are  sampling  points  of  the  2-transform. 

The  sampling  points  of  X(Z )  are  A,  B,  C,  D.  When  X(Z)  is  evaluated  at  those  points,  we  get 
the  DFT  X(K)  of  the  sequence. 

As  a  geometrical  interpretation,  we  see  that  the  exact  values  of  Z\,  Z2,  Z3  present  an  exact  "sym¬ 
metry"  towards  the  sampling  points  A,  B,  C,  D,  in  the  sense  that  the  product  of  the  magnitudes  of  the 
three  phasors  (from  the  sampling  point  to  the  zeroes  Zj,  Z2,  Z3)  gives  exactly  the  value  2,  for  each 
sampling  point.  In  Fig.  13,  we  sketched  the  three  phasors  for  the  sampling  point  B. 

For  point  B: 

10  -  Z,)  (J  -  Z2)  0  -  Zj)l  -  l«tl  •  M  •  M  -  2.  (75.1) 

For  point  A: 

|(1-Z,)(1-Z2)(1-Zj)|-2.  (75.2) 

and  similarly  for  points  C  and  D. 

This  property  (75)  is  evident  when  looking  at  the  Z  transform: 

X(Z)  -  1  +  Z~l  -  Z~2  +  r\  (76) 

and  substituting  directly  the  sampling  points  A,  B,  C,  D.  But  from  a  geometric  point  of  view,  it  is  a 
rare  combination  of  the  zeroes  of  X(Z),  that  present  such  a  "nice"  symmetry. 
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Notice,  however,  that  these  specific  zeroes  of  Barker  code  4  (Zt,  Z2,  Z3)  do  not  present  the 
above  symmetry  towards  any  number  of  uniformly  distributed  sampling  points  on  the  unity  circle.  For 
example,  for  eight  uniformly  distributed  sampling  points,  one  of  them  will  be  Z  -  eJw/*,  and  substitut¬ 
ing  it  in  (76): 

I X(Z  -  eJw/i)  |  -  |1  +  e~JwlA  -e~J*n  +  e~JiwlA\  *  2, 

which  means  that  these  specific  zeroes  of  Barker  code  4  cannot  be  'used"  for  generating  higher  length 
codes. 

Of  course,  the  same  analysis  holds  for  a  complex  sequence  U,|  •»  1,  except  that  the  N  —  1  zeroes 
of  X(Z)  will  not  be  in  conjugate  pairs.  But  again  for  a  G-polphode,  these  N  -  1  zeroes  of  X(Z)  are 
required  to  present  the  above  "symmetry"  towards  the  AT  sampling  points  Z  -  WK. 

One  might  suspect  that  some  uniform  distribution  of  the  zeroes  of  X(Z)  will  give  the  desired 
symmetry.  A  moment  of  reflection  shows  that  it  is  impossible  since  we  have  N  —  1  zeroes  of  X(Z) 
and  A  sampling  points. 

This  means  that  if  there  is  a  solution,  the  zeroes  of  X(Z)  will  be  distributed  on  the  Z  plane  in 
some  "rare"  combination  (and,  of  course,  not  on  the  unity  circle). 

Beyond  the  above  "symmetry"  these  N  -  1  zeroes  of  X(Z)  must  fulfill  other  requirements. 


i 

1 


! 


Zn-v 


Suppose  we  found  some  "symmetric"  structure  (in  the  above  sense)  of  the  zeroes,  Zt,  Z2,  .... 


Then: 

~  Zx)  (Z  -  Z2)  . . .  (Z  -  Zjv_,)  Xx  ZNJX  +  X2  z"-1  +  . . .  +  XN 
X\Z)  ■*  ~v_i  “*  (77) 


XU)  ZN~l 


ZAM 

-  (-Z,)  (-Z2)  ...  (-ZN„X)  -  XN 


ZN-\ 


(78) 


Z-0 


i.c.; 

|Zj  •  Z2  •  •  Z^-\\  “  1,  (78.1) 

which  means  that  some  of  the  zeroes  are  outside  the  unity  circle  while  the  others  are  inside,  such  that 
their  product  has  unity  magnitude. 

Another  point  to  mention  is  that  the  necessary  condition  is  "similar"  to  designing  an  exact  all  pass 
discrete  filter  whose  finite  impulse  response  is  h (a)  -  (xi,x2,  ....  *„},  where  |x,|  -  1. 

In  Ref.  S,  it  is  shown  that  an  all-pass  discrete  filter  has  a  Z  transform  that  factors  to  terms  of  the 

form: 

H,lz)  ~  (,,) 

where  0  <  a  <  1  (a  is  real),  such  that  the  pole  and  the  zero  in  (79)  give  a  constant  amplitude  for 
every  frequency.  This  is  actually  more  than  we  need,  since  our  requirement  is  to  have  constant  ampli¬ 
tude  /only  in  the  N sampling  points  of  X(Z).  But,  of  course,  in  our  problem,  since  we  have  a  finite 
length,  we  don’t  have  poles  of  X(Z)  (except  those  in  the  origin),  and  we  cannot  get  terms  of  the  form 
(79). 
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FINAL  COMMENTS 

1.  It  is  interesting  to  note  that  Frank  codes  of  even  length  meet  the  requirements  of  constant 
amplitude  DFT,  and  R(K)  4-  R*(N  -  K)  -  0,  but  still  they  don’t  form  G-polphodes.  For  example, 
the  Frank  code  of  length  N  -  16  is: 

<£,(deg):  0,  0,  0,  0,  1  0,  90.  180.  -90.  i  0,  180.  0.  180,  |  0,  -90,  180,  90  |  (80) 


1  /o 

1/0 

1  /0 

w 

\ 

1  /— 90 

1  (0 

-W 

1  io 

1  /90 

w 5 

1  /180 

1  /90 

1  t?o 

- 

w* 

1  JO 

1/0 

1  /180 

-  w 

1  io 

1  1^90 

1  /180 

W 

1  /o 

1 Z2  „ 

1  £ 90 

-W5 

i  /iso 

1  /90 

i  (% 

-  w5 

where  W  -  e~Jlw/li  -  1/-22.50. 

Also,  it  is  easy  to  verify  that  R(K)  +  R*(N  -  K)  -  0,  but  clearly  some  time  sidelobes  of  the 
Frank  code  are  bigger  than  unity  magnitude. 

This  provides  evidence  again  that  our  analysis  gave  necessary  conditions,  but  not  sufficient  ones. 
Therefore,  we  have  to  search  for  the  solution. 

Note  also  that  Barker  4  codes  are  actually  a  special  case  of  Frank  codes.  The  analysis  can  help  in 
searching  for  structures  of  either  the  code  sequence  xh  or  the  distribution  of  the  zeroes  of  X(Z). 

2.  An  issue  to  be  further  investigated:  Is  it  possible  to  approximate  the  requirement  of  constant 
amplitude  DFT  of  the  sequence,  and  thus  approach  the  "Barker  level"  of  the  autocorrelation  function? 
At  least  intuitively  we  might  think  that  a  constant  amplitude  DFT  is  a  "good  property." 

CONCLUSION 

The  motivation  for  the  analysis  was  to  find  a  finite  length  code  with  Barker  type  autocorrelation 
function  beyond  the  known  ones.  Though  no  specific  code  was  found,  the  analysis  derived  necessary 
requirements  for  even-length  Barker  codes  and  G-polphodes.  These  requirements  can  reduce  the 
search  problem  for  the  above  codes. 

The  different  points  of  view  presented  here  (time  domain,  DFT  of  sequences,  the  Z-transform 
geometrical  interpretation,  linear  algebra)  might  also  suggest  structures  and  properties  of  good  codes, 
which  only  approach  the  Barker  level  in  the  autocorrelation  function. 
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Appendix 

FURTHER  PROPERTIES  OF  AN  EIGENVECTOR  G-POLPHODE 

We  have  seen  that  an  eigenvector  G-polphode  must  fulfill  Eq.  (68.2)  of  the  text: 

BX~  ±  1  X  (Al) 


Let  Xa,  Xb  be  eigenvectors  of  B,  which  correspond  to  X  -  +1,  X  -  -1,  respectively. 

BXa  -  Xa  (also  B*Xa  -  X,)  (A.2.1) 

BXb  -  Xb  (also  B%  -  -Xb)  (A.2.2) 

From  (A.2.1),  multiplying  both  sides  by  B. 

BBX,  -  BXa+  B'%  -  X, .  (A.3) 

Similarly,  from  (A.2.2): 

BBXb  -  -BXb+  B1  Xb  -  Xb,  (A.4) 


i.e.,  X„  and  Xb  are  also  eigenvectors  of  B2,  both  correspond  to  the  eigenvalue  X  -  1  of  B2. 


The  matrix  B J  (N  x  Ama‘rix)  is: 

1  0  0  0  •  •  •  0  0 

0  0  0  0  ••  •  0  1 

0  0  0  0  •  •  •  10 

S2-  • 

0  0  0  1  •  •  •  0 

0  0  1  0  •  •  0 

0  1  0  0  ••  •  0 

The  matrix  B 2  has  N  eigenvalues;  some  of  them  are  X  -  1 .  and  the  others  urc  X 
way,  the  eigenvalues  X  -  ±  1  of  B  map  to  the  eigenvalue  X  -  1  of  B2). 


(A.5) 


-1  (by  the 


From  (A.3),  (A.S): 


(A.6) 


t 


NRL  REFORT  8586 


i.e.,: 

XN  *“  Xlt  ■*  X},  Xf/~ 2  “  X4,  etc,  (A.7.1) 

The  same  analysis  holds  for  Af4,  so  that  if  a  G-polphode  eigenvector  exists,  it  must  have  the  struc¬ 
ture: 


i/«t 

1  lh 

1  jot} 

1  lh 

1  {ai 

'lh 

1/04 

'lh 

1  /“W 2+1 

1  /0W/2+I 

1  /2! 

1  lh 

1/23 

1  lh 

1/22 

1  lh 

Now  from  (A.2.1)  (by  conjugating)  we  get: 

X •  -  X\,  (A.8) 

BB •  "►  4T.  - 1 *«•  (A.9) 

and  similarly, 

BX\--X\.  (A.  10) 

Equations  (A.2.1)  and  (A.9)  mean  that  if  Xt  is  an  eigenvector  of  B  (for  X  -  1),  then  Af*  is  also 
an  eigenvector  of  B  (also  for  X  -  1).  If  X,  is  real,  then  A'*#  -  X,  (they  are  identical). 

But  we  look  for  an  eigenvector  code  X,  (in  which  I  A’,  I  -  1).  We  have  seen  In  the  text,  that  such 
an  Xa  cannot  be  real  for  /  >  2.  So  Af„  if  it  exists,  is  a  polyphase  code.  Therefore,  is  another 
(linearly  independent  to  £,)  eigenvector  code. 

Similarly,  if  X>  is  an  eigenvector  code  (for  X  -  then  A'»  is  another  (linearly  independent) 
eigenvector  code. 

To  summarize,  the  eigenvector  G-polphodes  (for  /  >  2)  will  fulfill: 
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BXa 


BX> 


1  Zai 

'Hi 

1  /-Ofi 

1  /-a, 

1  [02 

i  Hi 

1  /-a2 

1 

f-ot  2 

1  /a. 

'Hi 

1  /~«3 

1  /-ay 

• 

,  BX  *  -  B 

• 

1  Hn/2+1 

1  H  A/2+1 

1  /-°W2+I 

i  /— aAf/2+i 

1  /a. 

*  Hi 

'till 

1  t2l 

'  Hi 

'Hi 

'tn 

1  l~a * 

'H i 

'Hi 

1/-0. 

'th 

'lh 

'Hi 

'th 

'th 

'Hi 

'Hi 

'  th 

'  th 

1  Hm* » 

m 

1  /Pm*i 

.  m  -  b 

l  /~fim*\ 

- 

1 

1  Hi 

'Hi 

'til 

'th 

'Hi 

'Hi 

» Hh 

'th 

(A.ll) 


(A.12) 


Note  that  in  (A.ll),  and  similarly  in  (A.  12),  out  of  the  S  equations  for  Jj,,  we  have  S/2  -  I 
redundant  equations  which  can  be  erased;  the  last  equation  is  identical  to  the  second,  the  (S  -  l)th 
equation  is  identical  to  the  third,  etc.  Unis,  if  we  erase  the  last  S/2  -  I  equations,  we  are  left  with 
S/2  +  I  equations  (some  of  them  are  also  redundant)  that  should  be  solved  parametrically  for  a,,  <*>. 

<»i . «*/)• 

Since  we  have  parametric  relations  for  the  <*,'s,  we  have  to  choose  some  of  them  such  that  we  get 
the  desired  autocorrelation  function.  This  again  is  a  search  problem. 

Applying  the  above  analysts  for  A’  -  16  we  get: 

BX.  -  X. 

XJ  -  II 

1/0,1. 
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After  solving  the  N/2  +  1-9  equations,  we  get  the  parametric  relations: 
(1  l a 3  +  1  [<Xi  +  1  /a?!  ™  1/2  (1  +  1 

11  [0l2  +  1  /f*4  +  1  /jH  +  1  /«g]  "  (1  [Oi  ~  1  /(*()], 

(1Z^~1Zai-1Z^+1/£«“(1Z^~1Za z1  1 


and 


cos  (2  '2ir/16)  ’ 


(1  (a 4  -  1  /o6l  -  |1  ^£2  ~  1  A»gl 


1  -  cos  (2rr/16) 
cos  (3  •  2w/16) 


(A.  13) 


Since  the  conjugate  code  X*„  is  also  an  eigenvector  (also  for  X  —  1),  we  must  have  the  above  four 
relations  when  1  (a^  is  replaced  by  1  /-o^.  Adding  and  subtracting  equations  in  the  above  eight  rela¬ 
tions  (complex)  gives  eight  real  relationslwith  cos  a,  and  sin  a,). 

For  example,  the  first  equation  of  (A.  1 3)  together  with  the  corresponding  one  (with  1  /-a,) 
result  in: 

cosaj  +  cos  a j  +  cos  a7  -  1/2  (cos  a)  +  cos  a*), 
and 

sin  <*1  +  sin  aj  -  sin  ai  —  1/2  (sin  <*|  +  sin  a*). 


A  code  (of  length  16)  that  satisfies  the  eight  relations  is  a  candidate  for  eigenvector  G-polphode 
(has  to  be  verified  in  the  time  domain). 

Note  also  that  if  the  structure  (A.U)  Is  a  G-polphodc.  it  is  required  that  for  any  length  M 

2tr/3  <  (a}  -  0|)  <  4tr/3f 

since  the  second  step  of  the  autocorrelation  process  gives: 

2)-  1/0+  1/oj-a,. 

and  its  magnitude  must  be  smaller  than  or  equal  to  one. 


